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Представлены результаты численных расчетов турбулентного закрученного осесимметричного течения. Расчет турбулентных ха- 
рактеристик проведен на основе дифференциальной модели турбулентности. Приведены результаты сравнения численных рас- 
четов с экспериментальными данными. Показаны результаты численного исследования о влиянии геометрии рабочей зоны 
пневматического классификатора на характер течения. 


Центробежные аппараты широко используют во 
многих отраслях народного хозяйства. Вопросы про- 
ектирования новых моделей и совершенствование су- 
ществующих конструкций тесно связаны с анализом 
процессов, происходящих в этих аппаратах. Гидроди- 


намическая картина течения, реализуемая в таких 
устройствах, оказывает существенное влияние на весь 
технологический процесс. Очень часто в пневматиче- 
ских аппаратах в качестве рабочей среды используют 
воздух, а сами конструкции имеют сложную геоме- 

109 


Известия Томского политехнического университета. 2005. Т. 308. № 3 


трическую форму. При этом режимы течения всегда 
являются турбулентными. Физическое моделирова- 
ние таких процессов связано с большими затратами 
средств и времени. Методы математического модели- 
рования можно рассматривать как один из перспек- 
тивных способов успешного решения задачи о нахож- 
дения поля скорости, при этом получение аналитиче- 
ских решений практически невозможно или связано с 
низкой степенью достоверности. Следовательно, 
единственным способом решения поставленных за- 
дач можно считать метод численного моделирования. 

В работе рассматривается центробежный аппа- 
рат, который представляет собой неподвижный ци- 
линдрический корпус с осевым выходным патруб- 
ком в верхней части. Закрученный турбулентный 
поток подается через кольцевое отверстие в ни- 
жнем основании рабочей камеры. 

Математическое моделирование турбулентного 
течения в таком аппарате основывается на осреднен- 
ных уравнениях Навье-Стокса, в которых рейнольд- 
совские напряжения записываются с использовани- 
ем градиентной модели Буссинеска. Эти уравнения в 
цилиндрической системе координат имеют вид [1]: 
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Представленная система уравнений Рейнольдса 
является незамкнутой. Для замыкания этих уравне- 
ний необходимо использовать модель турбулентно- 
сти. Несмотря на значительные успехи современ- 
ных исследователей в вопросах моделирования тур- 
булентности, включая прямое численное модели- 
рование и методы крупных вихрей, до сих пор наи- 
более оптимальной по возможностям использова- 
ния в инженерных расчетах считается двухпараме- 
трическая дифференциальная модель турбулентно- 


сти. В развитии такого подхода авторами была соз- 
дана оригинальная трехпараметрическая модель 
турбулентности на основе уравнений переноса для 
кинетической энергии турбулентных пульсаций, 
удельной скорости её диссипации и кажущейся тур- 
булентной вязкости [2]. Эта модель, записанная в 
цилиндрической системе координат, имеет вид: 
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Здесь константы модели турбулентности имеют 
следующие значения: 


С, =0,09; С 2 =1- 


О’о.Г . 


С 3 = 1 — у = ; ст, =ст ш =0,35; * = 0,4. (3) 


Математическая модель турбулентного течения 
решалась численно с использованием метода конеч- 
ных разностей. Решение строилось в переменных 
скорость - давление. При записи исходных уравне- 
ний применялась схема физического расщепления 
по времени [3]. Разностный аналог исходной систе- 
мы дифференциальных уравнений в частных произ- 
водных записывается на разнесенной сетке, что 
обеспечивает высокую точность получаемого реше- 
ния. Представление конвективных членов уравне- 
ний переноса с помощью экспоненциальной зави- 
симости [4] позволило получить второй порядок 
точности решения по координатам и устойчивость 
расчетной схемы. Численное решение определялось 
на основе неявного метода переменных направле- 
ний с использованием алгоритма прогонки. 
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Для рассматриваемой задачи использовались 
следующие граничные условия. Во входном сече- 
нии все переменные имеют постоянные значения. 
На выходе из аппарата ставились мягкие условия 
Неймана. На твердых непроницаемых поверхно- 
стях имеет место условие прилипания жидкости и 
вырождение характеристик турбулентности, кроме 
удельной скорости диссипации. В этом случае зна- 
чение удельной скорости диссипации кинетиче- 
ской энергии турбулентных пульсаций определяет- 
ся из баланса между молекулярной диффузией и 
деструкцией диссипации. В зависимости от ориен- 
тации границы можно получить выражения: 
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где индекс м> соответствует координатам стенки. 

Для проверки работоспособности созданной мо- 
дели (1-4) и оценке достоверности получаемых ре- 
зультатов были проведены тестовые расчеты. Иссле- 
довалось течение между параллельными неподвиж- 
ными дисками с направлением потока от периферии 
к оси. На рис. 1 показано распределение радиальной 
компоненты вектора скорости поперек канала на 
различных радиусах от входа. В качестве масштаба 
скорости выбрана щ, радиальная скорость на входе в 
канал. За масштаб длины взята высота междисково- 
го пространства Н. Символами обозначены экспери- 
ментальные данные [5]. Цифры соответствуют зна- 
чениям г/г 0 : 1 - 0,6, 2 - 0,4, 3 - 0,275, 4 - 0,185, где 
г - текущий радиус, г„ - радиус входа. 
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Рис. 1. Сопоставление численных и экспериментальных дан- 
ных при Ке=1269 



1 — • — • — 

1 

— • — • — 

• 

1* / 

г* * 

1 1 

1 ■ ■ 

. 2 . 

і ■ 1 

* / 1 

1 ■ 

V ♦ 

• 4 4 

3 


♦ / 

V 



■ ♦ ♦ 

♦ ♦ 


+ * 

3 

> 4 

















\ + + і 


Хорошее согласование расчетных и эксперимен- 
тальных данных позволяет сделать вывод о возмож- 
ности использования построенной математической 
модели для численного исследования закрученных 
турбулентных течений в центробежных аппаратах. 


Численное исследование гидродинамики турбу- 
лентного закрученного течения в рабочей камере 
центробежного аппарата представлено на рис. 2 и 3. 

На рис. 2 изображено характерное распределе- 
ние линий тока в рабочей камере, из которого вид- 
но, что вблизи внешней стенке центробежного ап- 
парата формируется интенсивное рециркуляцион- 
ное течение. Основное течение входящей жидко- 
сти наблюдается вдоль вертикальной оси аппарата. 
На рис. 3 представлены линии одинаковых значе- 
ний окружной компоненты вектора скорости. 
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Рис. 3. 



Распределение изолиний окружной компоненты ско- 
рости 
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Рис. 4. Изолинии окружной компоненты скорости при Н/Ку 
а) 0,5,6) 0,1 

Как видно из рисунка, имеет место тенденция к 
существенному увеличению значений окружной 
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скорости по мере приближения потока к выходно- 
му патрубку. 

Одним из параметров, влияющих на гидродина- 
мическую обстановку в рабочей камере центробежно- 
го аппарата, является высота конструкции. Результа- 
ты исследования влияния этого параметра на распре- 
деление окружной скорости представлены на рис. 4. 

Как видно из представленных иллюстраций, 
уменьшение высоты рабочей зоны не ведет к суще- 
ственным изменениям поля окружной компонен- 
ты скорости. Однако следует отметить, что с изме- 
нением геометрии существенно меняется напра- 
вление движения основного потока. Таким обра- 
зом, угол между аэродинамической силой и цен- 
тробежной силой может существенно изменится. 
Это явление может оказать значительное влияние 
на процесс разделения тонкодисперсных частиц в 
центробежном аппарате. 

Созданная математическая модель турбулентно- 
го течения в центробежном аппарате может быть ис- 
пользована для оптимизации режимных и геометри- 
ческих параметров существующих устройств, а также 
при создании новых перспективных конструкций. 
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